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Spectroscopy has an illustrious history delivering serendipitous discoveries and providing a strin¬ 
gent testbed for new physical predictions, including applications from trace materials detection, to 
understanding the atmospheres of stars and planets, and even constraining cosmological models. 
Reaching fundamental-noise limits permits optimal extraction of spectroscopic information from 
an absorption measurement. Here we demonstrate a quantum-limited spectrometer that delivers 
high-precision measurements of the absorption lineshape. These measurements yield a ten-fold 
improvement in the accuracy of the excited-state ( 6 P 1 / 2 ) hyperfine splitting in Cs, and reveals a 
breakdown in the well-known Voigt spectral profile. We develop a theoretical model that accounts 
for this breakdown, explaining the observations to within the shot-noise limit. Our model enables 
us to infer the thermal velocity-dispersion of the Cs vapour with an uncertainty of 35ppm within an 
hour. This allows us to determine a value for Boltzmann’s constant with a precision of 6 ppm, and 
an uncertainty of 71ppm. 


Spectroscopy is a vital tool for both fundamental and 
applied studies. It is critical to drive improvements in 
both precision and accuracy in order to gain maximum 
physical information about the quantum absorbers - even 
as they are perturbed by the probe. The immediate mo¬ 
tivation for this work has arisen out of a call to the scien¬ 
tific community to develop new techniques to re-measure 
Boltzmann’s constant, /cb, in preparation for a redefini¬ 
tion of the kelvin [1]; however, advances in absorption 
lineshape measurement and theory will find applications 
in accurate gas detection and monitoring [2], studies of 
planetary atmospheres [3], thermometry in tokamaks [4] 
and understanding distant astrophysical processes [5, 6]. 
The accurate measurement of the natural linewidth and 
transition frequencies, which can be directly related to 
the atomic lifetime, level structure and transition proba¬ 
bilities, is important for testing atomic physics [7-9]. 

Here we present measurements of a transmission line- 
shape of cesium (Cs) vapor with a quantum-limited 
transmission uncertainty of 2 ppm in a 1 second mea¬ 
surement, which is to our knowledge, a factor of 16 times 
superior to anything previously demonstrated [10]. This 
extreme precision allows us to directly detect subtle line- 
shape perturbations that have not been previously ob¬ 
served. This observation prompted the development of 
a theoretical model, which now allows us to discriminate 
between the internal atomic state dynamics and their ex¬ 
ternal motional degrees. Using the model we are able to 
estimate the velocity dispersion of the atoms with a preci¬ 
sion (standard error) of 53 ppm during a single line scan, 
taking ~ 30 seconds. This is consistent with the the sam¬ 
ple standard deviation (also 53 ppm) over multiple scans, 
demonstrating the excellent reproducibility of the spec¬ 
trometer. The measurement precision averages down to 


3.7 ppm after 200 scans. These values are ~ 7 times bet¬ 
ter than the previous best results [10-12], and also yield 
a 10-fold improvement in the uncertainty of the excited 
state hyperfine splitting in Cs [13, 14]. Our measure¬ 
ment of the homogenous broadening component of the 
line shape has a precision within a factor of two of the 
best ever measurement of natural linewidth. Modest im¬ 
provements in the probe laser performance would deliver 
a new result for the excited state lifetime of Cs in a sys¬ 
tem that is experimentally and theoretically much sim¬ 
pler than that typically used for lifetime studies [15, 16]. 

At thermal equilibrium, the velocity distribution of 
atomic absorbers in a vapour cell is related to the temper¬ 
ature through the Boltzmann distribution. This simple 
and fundamental relationship forms an excellent founda¬ 
tion for a type of primary thermometry known as Doppler 
broadening thermometry (DBT) [17, 18]. DBT differs 
from the current leading technique for primary thermom¬ 
etry, which measures the temperature-dependent speed 
of sound in a noble gas contained in a well-characterized 
acoustic resonator. Until recently, the best determina¬ 
tion of /cb was made in 1988 with a total uncertainty 
of 1.7 parts-per-million (ppm) [19]. Refinements to this 
technique over the last twenty years have reduced its un¬ 
certainty to the 1 ppm level [20, 21], improving the uncer¬ 
tainty in the CODATA-recommended value for /cb from 
1.8ppm (CODATA-02) to 0.91 ppm (CODATA-10). De¬ 
spite these superb measurements, it is important that dif¬ 
ferent techniques are employed to measure /cb to reduce 
the possibility of underestimating systematic uncertain¬ 
ties that may be inherent to any single technique [22]. 

The DBT approach is one of the most theoretically 
transparent of the new approaches to primary thermom¬ 
etry. The spectrum of a gas sample is measured with 
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FIG. 1. Schematic of a high-accuracy linear absorption spectrometer based on a probe laser (P) that is tightly locked at a 
tunable frequency difference away from a master laser (M), which is itself stabilised to an atomic transition in Cs. The frequency 
stability of the master laser is compared to a reference frequency provided by an optical frequency comb (C). The probe is 
spectrally filtered with an optical cavity (OC), and actively power-stabilised with an acousto-optic modulator (AOM). The 
probe is then divided with high-precision into two beams using a Glan-Taylor (GT) polarizing prism followed by a Wollaston 
beamsplitter (W). One beam passes through a Cs cell embedded inside a thermal and magnetic shield and is measured by 
photodiode A (PD A); the other beam is measured directly by photodiode B (PD B). Temperature and thermal gradients are 
measured independently with calibrated standard platinum resistance thermometers (SPRTs). 


high precision, and a theoretical model is fitted to the 
measured data. If the model includes all of the relevant 
physics then one can accurately extract the contribu¬ 
tion to spectral broadening arising solely from the ther¬ 
mal distribution of atomic/molecular speeds. Ammonia 
probed by a frequency-stabilized CO 2 laser at 1.34 ^m 
was the first thermometric substance employed in a DBT 
experiment [23]. Subsequently, an extended-cavity diode 
laser at 2 yum was used to probe a ro-vibrational tran¬ 
sition of CO 2 [24]. In these first experiments, the line- 
shape was assumed to be a Gaussian or a Voigt profile (a 
Gaussian convolved with a Lorentzian). Since 2007, with 
the ambitious goal of approaching 1 ppm accuracy, sub¬ 
stantial experimental and theoretical improvements have 
been made to DBT using ammonia [25, 26], oxygen [27], 
ethyne [28], and water [29]. However, one key challenge 
peculiar to molecular absorbers is the need to account 
for complex collisional effects on the line shape [26]. We 
avoid this by using a dilute atomic vapour [30, 31] with 
a strong dipole transition for which a tractable, micro¬ 
scopic theory has been developed [32, 33]. 

The spectrometer is pictured in Fig. 1. The probe laser 
is spectrally and spatially filtered using a combination 
of an optical cavity of moderate finesse {F « 305) and 
single-mode fibre. This reduces the spontaneous emission 
content of the probe beam from 1.6% to below 0.01%. It 
is then delivered into a vacuum chamber in which the 
vapour cell is mounted in a thermal and magnetic shield. 
The temperature of the shield can be controlled to a few 
millikelvin and gradients are suppressed to the same level. 
The light is split into two output signals using a combi¬ 
nation of a Glan-Taylor polarising prism and Wollaston 
beam-splitter. The ratio of powers in the output beams 
is stable to better than 10“®. One output directly illu¬ 
minates a photodiode to give the incident power, whilst 
the other passes through the vapour cell and is then de¬ 
tected. The ratio of these photodiode signals gives us 
the transmission ratio, T. The incident power is actively 


controlled while the frequency of the incident light is set 
to an absolute accuracy of 2 kHz. The temperature of 
the vapour cell is monitored directly using a calibrated 
capsule-type standard platinum resistance thermometer 
(CSPRT), allowing us to quantify the performance of our 
Doppler thermometer. Further technical details about 
the apparatus are given in the supplementary informa¬ 
tion (SI). 

We measure the transmission through the vapour cell 
as the probe laser frequency scans across the Gs D1 tran¬ 
sition ( 6 S 1/2 - 6 P 1 / 2 ), shown in Fig. 2(a). The relative 
noise in the measurement of the atomic transmission is 
just 2 ppm in a 1 s measurement (at the highest opti¬ 
cal powers used). The two dips seen on Fig. 2(a) come 
from the hyperfine splitting in the excited state and our 
excellent signal-to-noise ratio enables us to determine 
this splitting to be /hfs = 1167.716(3) MHz consistent 
with previous measurements (/hfs =1167.680(30) MHz 
[13, 14]), but 10 times more accurate. 

The absorption line-shape, 7It(/), of an atomic tran¬ 
sition with a rest-frame transition frequency /o, is 
given by a convolution of the natural (half) linewidth 
Fi = l/(47rr) due to the finite lifetime r of the tran¬ 
sition, with a Gaussian distribution having a e“^-width 
r£) = /o \/‘2kBT/{mcsC^) due to Doppler broadening of 
atoms of mass m at temperature T. By fitting a line- 
shape to the measured transmission data, with FD as 
a fitting parameter, we extract the Doppler component, 
from which we infer the temperature. Systematic er¬ 
rors in F^* are quantified by converting the indepen¬ 
dent PRT temperature measurements into an expected 
Doppler width, F™"^, using the CODATA values for fee 
and mcs, which we compare with FI). 

If the natural lineshape of the transition is a 
Lorentzian, then the resulting convolution is the well- 
known Voigt profile, which is commonly used to model 
dipole resonances. However, our extremely low-noise 
transmission measurements reveal deviations from the 
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FIG. 2. (a) Measured transmission at 296 K for the highest 
intensity, 1/Ias.t = 0.0028, averaged over 200 scans [points] 
and theory [solid], (b) Residuals including Voigt only [dark- 
green] and first-order correction [light-blue], (c) Etalon pre¬ 
factor, Eq. 1, with 6 etalons [solid] and data divided by 
second-order power-corrected Voigt proHles [points], assum¬ 
ing Fn = 2.327 MHz. (d) Residuals after including first-order 
correction with 0 etalons [light-blue] and 6 etalons [dark-blue], 
and second-order with 6 etalons [dotted-red]. 


Voigt profile. Some of these deviations are technical in 
origin e.g. instrumental broadening due to the lineshape 
of the probing laser, residual spontaneous emission from 
the probe laser, unwanted optical etalons, and photode¬ 
tector linearity; however, there is an important funda¬ 
mental effect resulting from frequency-dependent optical 
pumping, which perturbs the atomic natural line-shape 
away from a Lorentzian. All of these effects, whether 
technical or fundamental, cause systematic perturbations 
to the line shape, and must be taken into account to 
model the lineshape accurately. 

To demonstrate these deviations, we fit a model con¬ 
sisting of two Voigt profiles separated by /hfS) to the 
raw transmission data shown Fig. 2(a). The residuals of 
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FIG. 3. Fractional deviation between the fitted Doppler width 
and that inferred from PRT temperature measurements at 
296 K. Green points are fits using the power-independent 
Voigt-only profile. Blue points are fits using the first-order 
intensity-dependent correction. Red points are fits using the 
second-order intensity-dependent correction. Solid curves are 
fits using our central estimate of Fl = 2.327 MHz; dashed 
curves are for Fr, = 2.320 MHz (upper) and Fi, = 2.334 MHz 
(lower) representing F_l ± oTj^ : the shaded regions indicates 
the sensitivity of Ff) to Fl in this range. Error bars represent 
standard errors of fits to the average of 200 scans for fixed F l ■ 


this Voigt-only fit have characteristic M-shaped features 
near resonance, with amplitude ^ 200 ppm at the highest 
probe powers, as shown in the green trace in Fig. 2(b). 
These features demonstrate the breakdown of the Voigt 
profile, indicating additional, unaccounted physics, which 
causes a spurious linear dependence of F^ on probe in¬ 
tensity, as shown in the green curve on Fig. 3. In what 
follows, we describe how we sequentially include addi¬ 
tional physics in our transmission model to remove these, 
and other, systematic effects, leaving only the shot-noise 
limit of our detection apparatus. 

In earlier work we showed that the M-shaped features 
in Fig. 2(b) arise from atomic population dynamics in¬ 
duced by the probe laser, which are significant even at 
exceedingly low intensities [32]. We subsequently calcu¬ 
lated corrections to the Voigt profile up to linear-order 
in the probe intensity [33], to account for optical pump¬ 
ing. Including these first-order intensity-dependent cor¬ 
rection in our model suppresses the features substantially, 
as shown in the blue trace in Fig. 2(b). 

Far from the resonances, oscillations are evident in 
Fig. 2(b). These arise from low-finesse etalons in the op¬ 
tics used to deliver the light to the atoms. The dominant 
etalon has an amplitude around 40 ppm, corresponding to 
interference between surfaces with power reflectivities of 
< 10“®. Although obviously technical in nature, reduc¬ 
ing the size of etalons beyond this already fantastically 
low level is an experimental challenge. Instead, we in¬ 
clude etalons in our transmission model, so that the total 
transmission is given by T = TatTet, where %t includes 
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Source 

Mr(fcs) (ppm) 

Statistical (from Fig. 4) 

5.8 

Lorentz Width (r|,‘)“ 

65 

Lorentz Width 

190 

Laser Gaussian noise 

16 

Optical pumping 

15 

Etalons (misidentification) 

15 

Etalons (unresolved) 

3 

Spontaneous Emission 

3.6 

Temperature 

1.9 

Temp. Gradient 

1.2 

PD Linearity 

1 

Zeeman Splitting 

<1 

Total (fit Fl) 

71 

Total (independent T^) 

191 


TABLE I. Error budget for the determination of the Boltz¬ 
mann constant at 296 K. Optical pumping shifts are reported 
at I/Isat = 3 X 10“® assuming second-order corrections to the 
Voigt profile. “Estimated from fits to scan data with free LW. 
^’Estimated from published uncertainty in Cs lifetime, plus un¬ 
certainty in the independent estimate of the laser linewidth. 
Ur denotes fractional uncertainty. 



Scan Number 


FIG. 4. Deviation between rfj* (using second-order inten¬ 
sity correction) and for each scan at 296 K, with 

r_L = 2.327 MHz. Colours and vertical lines demark differ¬ 
ent incident intensities (lowest at left, corresponding to nor¬ 
malised intensities in Fig. 1). Dashed lines are the mean of 
the corresponding set of scans, error bars are the estimated 
parameter error (Tier) for each selected scan, and are con¬ 
sistent with sample standard deviation. Etalons are included 
in each fit, with parameters fixed by fits to the average of all 
scans within a power level, i.e. etalon parameters are not free 
to vary within a particular power level. 


n etalons and a slowly-varying quadratic background: 

^ (fy a + Pf + jf 

n"=i (1 + aj sin(27r///j -b 0^))' 

Etalons are of particular concern for DBT since they 
introduce systematic errors in E^. We note that it is 
only because of our extreme transmission sensitivity that 
these features are revealed, and thus give the opportunity 
to suppress this systematic error. 

We add etalons to the model until the residuals far 
from resonance are consistent with a white-noise back¬ 
ground. In practice, it is necessary to include up to n = 6 
etalons for high power scans, with the smallest resolved 
etalons having amplitude ~ 7 ppm. The solid curve in 
Fig. 2(c) shows Tet for the highest power data; the points 
show the raw transmission data divided by the fitted Tit, 
demonstrating that Tit is uncorrelated with the atomic 
transmission profile shown in Fig. 2(a). 

The dark blue trace in Fig. 2(d) includes first-order 
power-dependent corrections, and 6 etalons. Away from 
resonance, the residuals are then consistent with the 
noise-floor of our apparatus. However, around the deep¬ 
est resonance we observe ^10 ppm features. To elimi¬ 
nate these we include second-order intensity-dependent 
corrections to the Voigt profile (see SI): 

"Tat ^ (l + - l)Vj2)(A)/w(A) 

+ - l)Vj2)(A)/K(A))' 

+ - 1)fJ3)(A)/K(A)), (2) 


where v = Tjj/Tl, A is the detuning from the res¬ 
onance, and the generalised Voigt profile is given by 
(A) = 7r-3/2 dx (i + and the 

conventional Voigt profile is 14(A) = yj^^(A). For a 
cell of length L with linear absorption coefficient a, the 
optical depth on resonance is p = aL, oc T is the lin¬ 
ear intensity-dependent coefficient [32], and < 72,3 oc T^ are 
quadratic intensity-dependent coefficients. 

The red dotted trace in Fig. 2(d) shows the residuals 
after including second-order intensity corrections. The 
residuals are suppressed to ~ 2 ppm, consistent with the 
detection noise floor, giving confidence that the trans¬ 
mission model accounts for systematic effects. 

Figure 3 compares the fractional deviation between 
Tp and F™"'" as a function of incident intensity for the 
three possible models for the atomic line shape: Voigt- 
only (green), first-order (blue) and second-order (red) 
intensity-dependent corrections. Naturally, the true T d 
should be intensity independent. The simplest Voigt-only 
theory clearly exhibits a linear intensity-dependence, 
leading to ~ 750 ppm discrepancy even at intensities 
as low as I/Isat = 2.8 x 10“^ where Igat = 2.5mW/cm^ 
is the saturation intensity for the transition [13]. From 
[33], we expect the first-order correction to have a resid¬ 
ual quadratic dependence on intensity, which is consis¬ 
tent with the blue curve on Fig. 3. Finally, the second- 
order correction (red curve) is seen to suppress all in¬ 
tensity dependence to below the measurement precision. 
The simultaneous removal of all systematic features from 
the residuals together with the elimination of intensity- 
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dependence in the recovered gives a high degree of 
confidence that all relevant physics is properly included 
in the theoretical model. 

Thermometry: In this section we demonstrate the 
power of our spectroscopic technique by applying it to the 
problem of primary thermometry, and determine a value 
for Boltzmann’s constant using atomic spectroscopy. 

The atomic vapour is brought to thermal equilibrium 
at a chosen temperature by suspending it within a ther¬ 
mal isolator with an independent means for temperature 
measurement accurate to 0.55mK (1.9ppm). After sys¬ 
tematically removing the effects of optical pumping and 
etalons, the largest source of uncertainty in comes 
from uncertainty in Fr,. We obtain a value for F^ di¬ 
rectly by fitting it to the data, which gives F|* = 2.327(7) 
MHz. This is consistent with, but also more precise than, 
an independent estimate, F™‘^ = 2.331(19) MHz, given 
by the sum of the natural linewidth of the Cs 6 P 1/2 
level, F|) = 2.287(6) MHz [13], and the laser linewidth, 
F'f® = 0.044(18) MHz. From Fig. 3 we estimate that a 
1 kHz change in F^ leads to ~ 5 ppm change in F^, so 
that the 6.5 kHz uncertainty in F^ contributes 32.5 ppm 
uncertainty in F|,‘ [65 ppm in Ur^ks)]. We note that if 
the probe laser linewidth were decreased to the kilohertz 
level, then our measurement would yield a state-of-the- 
art value for the excited state lifetime of Cs. 

We now briefly describe additional sources of uncer¬ 
tainty, which are quantihed in Table I (see SI for a com¬ 
prehensive description). Statistical error arising from the 
scatter in Fig.4 contributes 2.9 ppm to Ur(T^) [5.8ppm 
to Ur (fee)]. This is also consistent with the estimated 
standard error for each point in Fig. 4, shown as er¬ 
ror bars. The probe laser has gaussian noise of width 
0.88(39) MHz, contributing 8 ppm to Mr(Ff)) [16 ppm to 
Uriks)]- Residual effects of optical pumping after second- 
order power corrections contribute 7.5 ppm to Ur(FD), 
[15 ppm to Ur(fcB)]. Misidentification of etalon parame¬ 
ters contributes 7.5 ppm to Ur(r|))) [15 ppm to Ur{kB)]- 
Possible unresolved etalons contribute 3 ppm to Ur(k-Q). 
Residual spontaneous emission from the laser, after filter¬ 
ing by the optical cavity, contributes 3.6ppm to Ur{k^). 

To determine a value for k-Q, we take a weighted 
mean of F^‘/F)],^"^ extracted from fits using all the 
second-order corrected data shown on Fig. 3 (red points). 
We find fce = 1.380545(98) x 10“^^ J/K, where the 71 
ppm uncertainty is calculated in Table I. This is 
consistent with the recommended CODATA value of 
1.380 648 8(13) x J/K. 

In conclusion, we have developed an atomic absorption 
spectrometer that operates with a transmission measure¬ 
ment precision of just 2 ppm. This has revealed several 
new phenomena. In combination with a theory that cor¬ 
rectly treats the interaction between light and an effusive 
vapour we are able to explain all of the observed effects 
to a level of precision never before demonstrated. The 
power of our technique is demonstrated by measuring 


the line-shape parameters of a gas at thermal equilib¬ 
rium. Our reproducibility is exactly consistent with the 
quantum-limits of measurement [34] which gives us con¬ 
fidence that we have captured all of the relevant physics. 

With our unprecentented sensitivity, we have measured 
the transmission spectrum of Cs at 895 nm at the shot- 
noise limit. From these measurements, we derived a 
value for Boltzmann’s constant with an uncertainty of 
71 ppm, which is consistent with the recommended CO¬ 
DATA value, and the hyperfine frequency splitting of the 
&P 1/2 level with an uncertainty ten times better than 
ever previously reported: 

fee = 1.380 545(98) x lO'^^ J/K 
/hfs = 1167.716(3) MHz. 

Our current performance permits a measurement of the 
Boltzmann constant with a precision of 6 ppm after only 
a few hours of data acquisition. Our total uncertainty is 
dominated by the imprecision in the literature value of 
the excited state lifetime of the Cs D1 transition (known 
to 0.26% [13]). A future experiment designed to fur¬ 
ther reduce the amplitude of etalons, and using a more 
broadly tuneable (~ 30 GHz) and higher power probe 
laser with sub-kHz linewidth would eliminate probe laser 
effects, allow better estimation of etalon contamination, 
and lower the contribution of optical pumping. This 
opens the door to the best measurement of the Cs D1 
transition lifetime, and a ppm-level measurement of fce 
thus making an important contribution to efforts to re¬ 
define the kelvin. 
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APPARATUS 

Light Source: A dual-beam linear absorption spec¬ 
trometer operating at 895 nm was built using a pair 
of extended cavity diode lasers (ECDLs). The first of 
these, the master laser, is locked to the 6 S'i /2 F = 4 
to 6 Pi /2 F = 3 transition of the D1 line in Cs at 
894.6054 nm. This provided an optical reference fre¬ 
quency with a square-root-Allan-variance (SRAV) of 
~2kHz for timescales between Is and 30s, as shown by 
the grey circles on Fig. 1. These are the relevant time- 
scales for the spectroscopy, reflecting the time to take a 
single scan. The other laser is frequency locked to the 
master at a user-selectable frequency offset by stabilising 
the heterodyne beat frequency against a tunable radio 
frequency (RF) reference oscillator. This frequency lock 
contributed no additional frequency instability to this 
slave laser (red triangles). The resulting slave laser fre¬ 
quency stability is shown by the olive trace on Fig. 1. The 
RF reference oscillator was step tuned under computer 
control to ensure that the transmission data was syn¬ 
chronous with the selected frequency. The wavelength of 
the slave laser was centered on the twin absorption peaks 
at 894.5793 nm, corresponding to the 6 S'i /2 F = 3 to 
6 Fi /2 F = 3 and F = 4 transitions. The offset frequency 
was tuned in increments of 50 MHz over a 6500 MHz span 
that captured the atomic absorption features. 

The slave laser output was coupled into a Fabry-Perot 
resonator whose length was actively controlled to keep it 
resonant with the slave laser light. The output of the cav¬ 
ity was coupled into an acousto-optic modulator (AOM) 
which was used to actively stabilise the power in the de¬ 
flected beam. The output of the AOM was coupled into a 
single-mode fibre for spatial hltering and for delivery into 
the vacuum chamber. This fibre contained an auxiliary 
detector on a fibre coupler which provided the signal to 
stabilise the input power into the chamber. The rest of 
the light in the fibre was sent to a Cs sample cell embed¬ 
ded in a thermal and magnetic shield. 

Spontaneous Emission: In a separate experiment 
the slave laser light was tuned to the centre of the atomic 
resonance and passed through a heated Cs cell for which 
the optical depth was 30 i.e. Tat ~ ~ 10“^^. A 

measurement of the residual transmission demonstrated 
that the laser output contained 1 . 6 % broadband emission 
at its operating point. 

The laser light was filtered using the fundamental 
mode of a scannable Fabry-Perot cavity (F =305 and 
free-spectral range 17.6GHz). This was necessary be¬ 
cause broadband emissions contaminate the measure¬ 
ment by creating a transmission offset associated with 
the non-resonant photons. Low level feedthrough of light 



FIG. 1. Square-root-Allan-variance (SRAV) of the absolute 
frequency stabilities of the master laser (grey circles), slave 
laser (olive triangles) and the heterodyne beat (red triangles). 


from higher-order spatial modes of the optical cavity (~ 
20 % of the main mode) was heavily attenuated by cou¬ 
pling the beam into a single-mode optical fibre prior to 
introduction into the evacuated chamber. 

The combined cavity and single-mode filtering pro¬ 
vides a reduction in the spontaneous emission by a factor 
2 F/ 7 r « 195. We numerically modeled the effect of the 
residual spontaneous emission (0.008%) on the measured 
Cs spectrum and found that the perturbation to the fit¬ 
ted atomic Doppler width vjj was only 1.8 ppm (equiv¬ 
alent to 3.6ppm in ^^(/cb)). This systematic error can 
easily be reduced to below the 1 ppm level by a cavity 
finesse of over 1000 . 

Instrumental linewidths: The probe laser spectral 
lineshape can be expressed to a high degree of accuracy 
as a Voigt function [1] with a Lorentzian half-width-at- 
half maximum of P'f^® and Gaussian e“^-halfwidth of P*^®. 
When this instrumental function is convolved with the 
atomic absorption profile (which, to first order, is also 
a Voigt function with parameters Pi and Pd), the ap¬ 
parently observed atomic lineshape is also a Voigt func¬ 
tion with Lorentzian and Gaussian components given by 

Pi = and Pb, = y^(Pg )2 -b (P’g")^, respectively. 

The power spectral density (PSD), S^{f), of the probe 
laser frequency noise was directly measured by using the 
frequency-dependent absorption on the side of an atomic 
transition. Two regions of the PSD can then be identified 
[1]; at higher Fourier frequencies we see a dominant white 
frequency process, while at lower frequencies we observe 
steeper noise which in our case has a dominantly flicker 
(1//) character. White noise components (with ampli¬ 
tude So) produce a Lorentzian lineshape with HWHM 
given by ttSq /2. The integrated probe laser noise at lower 
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frequencies give rise to a Gaussian lineshape which has 
an estimated e-fold width of = 0.88(39) MHz when 
integrated over the 40 ms observation time for each fre¬ 
quency point in the scan. 

Optical power control: The spectrally purified light 
is passed through an acousto-optic modulator (AOM) 
and the deflected beam is coupled into a non-polarising 
fiber splitter, which splits the power in the ratio 9:1. The 
lower power beam is measured on a photodiode and a 
feedback loop adjusts the AOM input to maintain a con¬ 
stant optical power in the fibre. The higher power out¬ 
put is out-coupled and used for spectroscopy. The po¬ 
larisation of this beam is fixed by passing it through a 
Gian-Taylor prism, and then a Wollaston beamsplitter is 
used to divide it into a reference and sample beam. The 
sample beam passes through the Cs vapour and common¬ 
mode amplitude variations in this signal are removed by 
dividing it by the reference beam signal. 

Conventional thermometry: The Brewster-angled 
Gs sample cell is embedded in a cylinder of copper in¬ 
side a multilayered thermostat, which is itself inside an 
evacuated chamber. The temperature of the thermo¬ 
stat can tuned using a thermo-electric cooler attached 
to the outer-most layer of shielding. The temperature 
at the cell was monitored using a capsule-type standard 
platinum resistance thermometer (GSPRT) that was cal¬ 
ibrated to the ITS-90 temperature scale with a temper¬ 
ature precision of better than 0.5 ppm in the range be¬ 
tween 273 K and 300 K. However, our resistance meter 
(HP Model 3458A) limited our temperature measure¬ 
ment uncertainty to 1.2 ppm. Measurements using a sec¬ 
ond calibrated GSPRT located in the copper cylinder at 
a distance ~ 50 mm away from the primary GSPRT re¬ 
vealed no temperature gradients above the 1.2 ppm level. 

The ITS-90 temperature scale, Tgo, itself is known to 
deviate from the true thermodynamic temperature, T by 
T -Tqo = 3.2±0.4mK at Tgo = 296 K[2]. When the 
uncertainty in this correction is included, the total un¬ 
certainty due to conventional thermometry is 1.9 ppm. 

Magnetic Shielding: The magnetic environment was 
nulled around the Gs cell using a dual-layer shield made 
from high permeability metal. Numerical modeling and 
experimental characterisation of this shield demonstrated 
that the apparent broadening of the atomic absorption 
line due to Zeeman splitting of the otherwise degenerate 
levels contributed less than 1 ppm uncertainty in ks ■ 

Optical detection: Reverse-biased silicon photodi¬ 
odes were used to measure the optical power in each of 
the reference and sample beams. The photosignals were 
measured by digital multimeters (DMMs) recording the 
voltage generated by the photocurrent passing through a 
load resistor. The noise contributions to each measure¬ 
ment of the atomic transmission from shot noise, Johnson 
noise and the DMMs are shown in Fig 2 as a function 
of input signal level to each photodiode (the stair-case 
structure revealed in the figure arises from a change in 


the DMM range). It is clear that the transmission mea¬ 
surement are shot-noise limited in several bands of pow¬ 
ers that lie just near the maximum of each range. At the 
highest powers used in the experiments (^0.6/rW), shot 
noise contributed more than 50% of the total measure¬ 
ment noise. 

The linearity of the photodetection chain was mea¬ 
sured in two ways. In the first approach, the ratio of 
the transmitted and incident detector signals were mea¬ 
sured over the full range of common-mode input power 
variation. This provided a reliable calibration curve for 
which the deviations from linearity at high powers could 
then be corrected. For low powers, this method did not 
provide sufficient signal-to-noise for calibration. Instead, 
we monitored the ratio of absorption coefficients for the 
two optical transitions. The second approach relies on 
the fact that the Glebsch-Gordon coefficients for these 
transitions require that the linear (i.e. low-power) ratio 
of absorption coefficients should be 1:3 [3]. Forming this 
ratio, and enforcing that the extrapolation to zero power 
agrees with the 1:3 condition, we have a guaranteed, in 
situ and independent way to determine any nonlinear¬ 
ity in the detector response at low powers. In this way, 
we were able to resolve and remove the nonlinearities in 
the photodetection scheme. The technique allowed us to 
measure offsets corresponding to < 8 pW allowing us to 
ensure linearity over a dynamic range of 10^ : 1. A more 
complete description of this method is found in Ref. [4]. 

ERROR BUDGET 

Here we give a detailed description and derivation of 
the uncertainties reported in Table I of the main text. 
Ur{X) denotes the fractional uncertainty in X. 

Statistical: For each power level, the statistical con¬ 
tribution is computed by fixing the etalon parameters 
(to values obtained from fits to coaveraged scans, see be¬ 
low), and extracting from each scan independently. 
These are shown in Fig. 4, and we find a statistical error 
of 2.9ppm in Ff), [which represents 5.8 ppm in Uriks)]- 

Lorentz Width: Our probe laser has a Lorentzian 
linewidth, which adds to the atomic linewidth, F^ = 
Ff) -I- F'))®. To determine F l, there are two alternatives: 
either to treat F^ as a free fitting parameter (together 
with F^)), or to fix F|) to the best reported literature 
value and independently estimate F”^®. 

In the former route we fit to the average of 200 in¬ 
dependent scans for each probe power level we retrieve 
a value of F|* = 2.327(7) MHz. In the latter route, 
we take F|) = 2.287(6) MHz for the natural linewidth 
of the Cs 6 Pi /2 level [3], and an independent estimate 
of the laser noise F”))® = 0.044(18) MHz (see above), 
which totals F“'^ = 2.331(19) MHz. These two estimates 
for Fi are consistent within their measurement uncer¬ 
tainty, and so we take the value with smaller uncertainty 
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FIG. 2. Top: Calculated shot noise (S) and Johnson noise (J) 
voltages as a function of the photosignal. The DMM (D) input 
noise was empirically determined and found to depend on the 
measurement range setting (indicated by the shaded regions). 
Bottom: Measured noise in the atomic transmission channel 
(circles) compared to the calculated (with no free parameters) 
contributions from the DMM, Johnson noise, shot noise and 
total (T) are also shown for comparison. 


Fl = 2.327(7) MHz. In particular, the solid curves in 
Fig. 3 assume the central value F^ = 2.327 MHz, and the 
dashed curves correspond to the upper and lower bounds 
of the uncerainty range, F^ = 2.320 MHz and 2.334 MHz. 

We note that the relative precision in our measurement 
of F^, is within a factor of two of the precision of the best 
measurements of the Cs D1 transition lifetime, rcs [5, 6], 
and matches the precision of the best synthesised value 
of Tcs [3] formed by combining these measurements. 

We also see from Fig. 3 that a 1 kHz variation in 
Fi leads to ^ 5 ppm change in F^, which is consis¬ 
tent with a 1:1 tradeoff between increases in F/, and de¬ 
creases in Ff)*, when Fd « 215 MHz. This implies that 
our ±6.5 kHz uncertainty in F®* contributes 32.5 ppm to 
Ur(Ffj‘) (equivalent to 65 ppm to Ur{kB))- The ±19kHz 
uncertainty in F“'^ represents 190 ppm in Ur{kB). 

Laser Gaussian Noise: The probe-laser spectrum 
has a Gaussian contribution [1]. This was measured and 
found to have a width of 0.88(39) MHz. The Gaussian 



Normalised Intensity (///sat) 


FIG. 3. Fractional deviation between the fitted Doppler width 
and that inferred from PRT temperature measurements at 
296 K. Green points are fits using the power-independent 
Voigt-only profile. Blue points are fits using the first-order 
intensity-dependent correction. Red points are fits using the 
second-order intensity-dependent correction. Solid curves are 
fits using our central estimate of Fl = 2.327 MHz; dashed 
curves are for Fr, = 2.320 MHz (upper) and Fi, = 2.334 MHz 
(lower) representing F_l ± oTj^ : the shaded regions indicates 
the sensitivity of rf)‘ to Fl in this range. Error bars represent 
standard errors of fits to the average of 200 scans for hxed F l ■ 


component sums in quadrature with Fd, and so con¬ 
tributes 8 ppm to Ur(F|,*) [16ppm to Ur{kB)]. 

Optical Pumping: The second-order power- 

dependent corrections suppress the fit residuals in (shown 
in Fig. 2(d) of the main text) to the measurement noise 
floor, so there is no statistically significant signature of 
optical pumping. Nevertheless, any residual intensity de¬ 
pendence in the second-order intensity-corrected Ff) is 
expected to scale as (7/Fsat)^- We therefore estimate the 
error arising from uncorrected third-order optical pump¬ 
ing by fitting a cubic through the highest intensity red 
data points in Fig. 3. This gives a maximum deviation of 
15 ppm in itr.(FD^) at the maximum probe intensity, and 
smaller contributions at lower intensities. Forming the 
weighted average over all intensities in Fig. 3 contributes 
7.5 ppm to Mr(F^*), [15 ppm in Mj.(^b)]- 

Etalons: Etalons enter into the uncertainty budget in 
two ways. Firstly, etalon parameters may be misiden- 
tified during fitting, and so contaminate F^. Secondly 
unresolved etalons masked by measurement noise may 
introduce residual systematic shifts to F^j. 

In Fig. 3 we fit Ff)^ with free etalon parameters. The 
weighted average of the standard errors is 8 ppm in 
Mr(Ff)‘), which is a convolution of the statistical noise 
discussed above and uncertainties arising from possible 
misidentification of the etalon parameters during fits. 
When the statistical errors in F^* are deconvolved from 
the standard errors in Fig. 3, we find the systematic error 
due to free etalon parameters in the fitting routine to be 
(8^ — 2.9^)^/^ = 7.5 ppm in ^^(Ff)*), [15 ppm in Ur{kB)]- 
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FIG. 4. Deviation between rfj^ (using second-order Voigt 
correction) and for each scan at 296 K, with 

rL = 2.327 MHz. Colours and vertical lines demark differ¬ 
ent incident intensities (lowest at left, corresponding to nor¬ 
malised intensities in Fig. 1). Dashed lines are the mean of 
the corresponding set of scans, error bars are the estimated 
parameter error (iln) for each selected scan, and are consis¬ 
tent with sample standard deviation. 


We bound the error arising from unresolved etalons in 
two ways: by synthesising data with false etalons, and by 
computing the fractional shift in after the addition 
of the etalon in the fit function. Both approaches 
give worst-case shifts in that are comparable to the 
amplitude of the etalon. The largest unresolved etalons 
are smaller than the measurement noise floor, so these 
cannot contribute more than 3 ppm to Ur{kB)- 


FIG. 5. Generalised Voigt profiles versus detuning, u — 100. 


the Rabi frequency experienced by the atom as it moves 
through regions of different intensity. We write the Rabi 
frequency in this form since in what follows we will ex¬ 
pand the atomic response in powers of fl^/(l + A^). For 
simplicity of presentation, we will assume that the equi¬ 
librium atomic level populations are Peq = {1/2,1/2, 0}, 
i.e. we ignore degeneracies in the ground state manifold. 
Accounting for degeneracies is a simple modification to 
the relative fractions for Pi and P 2 at equilibrium. 

The light is absorbed as it propagates through the 
atomic medium. While Eq. 1 is defined with respect to 
the rest-frame of a particular atom, the axial evolution 
of the field in the lab-frame is governed by [7] 


LINESHAPE CORRECTIONS AT SECOND 
ORDER IN INTENSITY 


Here we derive Eq. 2 of the main text, which gives the 
second-order intensity-dependent correction to the Voigt 
profile. We follow closely the derivation in [7]. 

We begin by computing the spectral dependence on the 
atomic populations for a three level atom, consisting of 
two ground states: one optically active, |1), and the other 
optically inactive, |2), and an excited state |3). Transi¬ 
tions between states |1) and |3) are optically driven, and 
state 13) can relax to either of the ground states. The 
population rate equations are [7] 


P = MP, 


M = 


1-I-A2 

0 

\ 1-I-A2 


2/3+ 

2(1-/3) 
nV(t)/2 

H-A2 


0 
0 

0 - 2 - 


( 1 ) 


where P = {Pi,P 2 ,P 3 }, A is the detuning between 
the laser frequency and the atomic transition, j3 is the 
branching ratio from state |3) to |1), is the peak atomic 
Rabi frequency (proportional to the electric field am¬ 
plitude) and \/fit) represents the temporal profile of 


d ln02 
dz 


2'k^I'^K / dA„ 


g-CA-A.jV^'" 

3/2/1 I A2 '1) 

7r+^(l-f A^J 


( 2 ) 

where P = P^ — Pi, K = po/J'^/Cv^^o^l'Co)) A* is the tran¬ 
sition dipole moment and uq = ^/ffkj^Pjrn is the mean 
thermal atomic velocity for an atom of mass m. The 
notation A„^ indicates that a given atom is subject to 
a detuning which depends on its axial velocity via the 
Doppler shift. If the field intensity is weak, so that on res¬ 
onance populations are negligibly perturbed from ther¬ 
mal equilibrium, then P = —1/2, and the integral yields 
D2(z) = ^ D 2 ( 0 )/Dg, we 

recover T = e~'^ KzV,,{A)^ which is the conventional 

Beer’s law result for a Voigt atomic lineshape. 


For larger perturbations to P become significant, 
and we need to compute corrections to P arising from 
the dynamics described by Eq. 1. Since we are concerned 
with small but significant pump induced perturbations to 
the thermal population, we do a perturbative analysis of 
Eq. 1 around the limiting case where D = 0. 


By inspection of Eq. 1 we see that only ever appears 
as part of the ratio D^/(l -f A^), so it must appear in the 
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same ratio in P, which expand in powers of [7]: 


v = Td/Ti^, and we define a generalised Voigt profile 


P(A) =00 + 01 


1 + A2 


+ 02 


1+A2 


2 

+ 0(176), (3) 


Vj")(A) 



g-(A-A„jVO 

7r3/2(l _|_ A2J" 


so that the conventional Voigt profile is 


where oq = — 1/2 and other aj are expansion coefficients, 
which contain all the complexity of the ballistic atomic 
motion across the beam, including averages over trans¬ 
verse velocity classes and trajectories across the beam. 
Here P is defined in the lab-frame, so implicitly depends 
on the axial and radial position (z,r). 

Substituting Eq. 3 into the LHS of Eq. 2, and evaluat¬ 
ing the integral over A„^ yields 


d In 172 
dz 


2 

2 (A) {e flf + 0(176), ( 4 ^ 

j=o 


where we have introduced a pseudo-parameter e = 1 
which we will use later to find a power-series solution. 


K(A) = Vji)(A) = Re{e-(^+*)'/"'erfc[(l - zA)/^]/v^}. 
( 2 ) 

V/ ^ may be evaluated by noting that 

1 r , 1 • 

(1 + a; 2)2 5 ™ \ (^2 _ 1 ) (^2^2 + !) + (!_ ^2) (^2 + 1) ^ 


Convolving this sum of two Lorentzians with a Gaussian 
thus yields a sum of Voigt functions, 


V/ 2 )(A) = lim 




I (e^ - 1) 

(A + i) 




rK(A) 


= Re{- 


' + 


a-e) 

(zA +1.2/2 _i)erfc[l^] 


TTl/'^ 


}■ 


Higher order correction terms can be calculated itera¬ 
tively, using the same method. In particular 


V/ 6 ) (A) = Re I (2iy (21A + 3- 2) + V 7 re-(^+*)"/'^"erfc [(1 - f A) /i^] (2(-2 + 3i A)i^^ - 4(A + if + /( 87 r:/ 4 )| . 


These are shown in Fig. 5. 

Equation (4) can be solved analytically for 172 ( 2 ;) (for 
any order of expansion), however it becomes cumber¬ 
some. Instead, we expand in powers of e, by substitut¬ 
ing 172(z,r) = 17 q( 2:,r) + £ 17f(z, r) + £2l7|(z,r) + ... into 
Eq. 4, and solving order-by-order in powers of e, with 
initial conditions 172(0,r) = 17Q(r) (where the radial de¬ 
pendence is parametric). Physically, this is simply an 
expansion around the zero-power limit. We find 


17^(z,r) 

17?(z,r) 

nl{z,r) 


nl{r)e 


-pvpA) 


n^{r) 


(5) 

, ( 6 ) 


176 (r) 


Q^g-pK(A)(g-pV.(A) _ l)v/ 2 )^^^ 

aoV;(A) 

'„2g-pV+A)(g-pV+A) _ 1)2-^V2)^^^2 




Q2g-pV+A)(g-2py+A) _ l)vd3)(A) 
2aoK(A) 


,(7) 


This gives Eq. (2) of the main text: 

T = e-pvP^) (l + q^i^e-pV-'A) _ 1)v/2)(a)/h,(A) 

+ - l)VYiA)/VfA)f 

+ - I)VJ6)(a)/h,(A)), 

where qi = (ail7Q)/(aol7Q), q 2 = (a2f7(i)/(agl7Q) and 
<73 = (a 2 l 7 (i)/( 2 aol 7 ^), and fij} = rdrn^(r). 

It follows that <71 is proportional to the intensity, whilst 
<72 and (73 are proportional to the square of the inten¬ 
sity. This is made explicit if we assume that the input 
probe beam has a gaussian transverse profile at z = 0 , 
172 (r) = 417ge“2(’’/’'o)^^ — 2ail7g/ao, 

q 2 = 16a2l7g/(3ag) and qs = 8 a 2 l 7 g/( 3 ag). In particu¬ 
lar, (72 = 4^2/3, so that (72 is determined by qi. Since we 
do not know the parameters ai_ 2 , or the optical depth we 
treat p and qi as fitting parameters for each resonance. 


where p = ogz > 0 . 

Since the intensity is proportional to 172(z, r), we com¬ 
pute the transmission z = L as 

^ /g°° rdrnfL,r) 

/g“rdrl72(r) ■ 
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